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Abstract. We study a higher order parabolic partial differential equation that arises in the context 
of condensed matter physics. It is a fourth order semilinear equation whose nonlinearity is the deter¬ 
minant of the Hessian matrix of the solution. We consider this model in a bounded domain of the real 
plane and study its stationary solutions both when the geometry of this domain is arbitrary and when 
it is the unit ball and the solution is radially symmetric. We also consider the initial-boundary value 
problem for the full parabolic equation. We summarize our results on existence of solutions in these 
cases and propose an open problem related to the existence of self-similar solutions. 
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1. Introduction 

We are interested in the initial-boundary problem for u = u{x,y,t) solving the 
following parabolic equation 

(1.1) dtu + = dei{D‘^u) + Xh, 

subject to the initial condition u{x,y.,0) = uo(x, y) and where h is some function 
depending in general on both space and time coordinates and belonging to some 
suitable Lebesgue space, A € M. This equation is to be solved for {x,y) G U C 
where U is open, bounded and provided with a smooth boundary, and t > 0. We 
will consider two different sets of boundary conditions: Dirichlet and Navier. 

The research has been supported by projects MTM2010-18128, RYC-2011-09025 and SEV-2011-0087, 
MINECO, Spain, and project RVO 67985840, Czech Republic. 
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This partial differential equation is a model that arises in the coarse-grained de¬ 
scription of epitaxial growth processes in the field of condensed matter physics. The 
function u = u{x, y, t) models the height of the growing film at the spatial point 
(cc, y) at time t. Note that the evolution is dictated by the competition among the 
determinant of the Hessian matrix of the solution and the bilaplacian. These terms 
model the dynamics at the solid surface. The function h models the introduction of 
new mass on the surface, and the parameter A measures the intensity of this income 
of new mass. A formal derivation of this model in terms of geometric quantities can 
be found in 0 and references therein. 

We have partially analyzed this model in a series of recent works |[T112] [3l ID . We 
have shown the existence of stationary solutions to this partial differential equation, 
in the case h were time independent, for two different sets of boundary conditions: 
Dirichlet and Navier. The proofs are different for the two sets of boundary conditions 
because the variational structure that is present in the Dirichlet problem is absent 
in the Navier one 0. For solutions that are radially symmetric we recover the 
variational structure in both cases, and the proofs of existence of solutions are built 
making an explicit use of this fact 0. Furthermore, in the case of radially symmetric 
solutions it is possible to prove non-existence of solutions for large enough data 0. 
For the evolution problem it is possible to prove local existence of solutions for 
arbitrary data and global existence of solutions for small (but otherwise arbitrary) 
data. Depending on the boundary conditions and the concomitant presence of a 
variational structure in the equation as well as on the size of the data it is possible to 
prove blow-up of the solution in finite time and convergence to a stationary solution 
in the long time limit 0. 

A summary of these results will be exposed in the next section. Together with 
these proven facts, there is a number of questions that remain open for both partial 
differential equation (fTTl i and its stationary counterpart. One of them is the existence 
of self-similar solutions and its possible role in the blow-up structure. We describe 
an open question related to the existence of self-similar solutions in section [3) 

2. Summary of previous results concerning existence of solutions 

In this section we list without proof some results we have recently obtained. The 
proofs can be found in the references quoted in the text next to the statement in 
question. We are going to consider the two following sets of boundary conditions u = 
dnU = 0 on O, which we refer to as Dirichlet boundary conditions, and u = Au = 0 
on which we refer to as Navier boundary conditions. Following 0 we find that 
the Dirichlet problem for the stationary version of (11.11) 

(2.1) A^u = det{D^u) + Xh, 
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where of course h is assumed to be time independent, has a variational structure: 


Theorem 2.1. The critical points of the functional 


Jx : Wg'^in) ->■ R 
(2.2) u Jx{u) 


— ( \/S.u'\^ dx dy — I UxUyUxy dx dy — X j hudxdy, 

2 Jn Jn Jq 


are weak solutions to the Dirichlet problem for (12. 11 1. 


This variational structure can be used to prove the existence of at least two solu¬ 
tions to the corresponding boundary value problem. 

Theorem 2.2. Let h € Then there exists a Aq > 0 such that for 0 < A < Aq, the 

Dirichet problem for equation (12.11) has at least two solutions in 

The proof can be found in H and makes use of the mountain pass geometry of 
functional dOi . For the Navier problem for equation (ED we cannot use the above 
mentioned variational methods because we do not know of any suitable functional 
in this case. Instead we have the following result: 

Theorem 2.3. Let h € Then there exists a Aq > 0 such that for 0 < A < Aq, the 

Navier problem for equation (12.11) has at least one solution in {Li) D (fl). 

The proof makes use of Banach fixed point theorem, see n. Note also that this 
proof can be immediately adapted for the case of Dirichlet boundary conditions. 

The radial problem corresponding to (12.11) reads 




r 

- [ru ) 
r 


- u'u" + Xh{r), 
r 


where r = x/x^~+^. In this case, for both sets of boundary conditions, the problem 
admits a variational formulation. And in both cases, due to the mountain pass 
geometry of the associated functional, it is possible to prove the existence of at least 
two solntions for small enough A, see JD- Furthermore, it is also possible to prove the 
non-existence of solutions for large enough A as well as rigorous bounds for the values 
of A that separate existence from non-existence, see 0. Note that these bounds are 
rather precise in certain cases when compared to the numerical estimations of the 
critical values of A calculated in El. With respect to the full evolution problem (ED, 
we can prove the following theorem: 

Theorem 2.4. Let T > 0; for any uq € IFg’^(D), any h G L^(0,T;L^(D)) and any 
A € R the Dirichlet problem has a unique solution in 

Xt ■■= C(0, T; IFg’^(D)) n L^(0, T; W^’^{n)) n W^’^{0, T; L^(D)), 
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provided T is sufficiently small. Furthermore, for any T € (0,oo] there exists a unique 
solution to this problem in the same space provided ||uo||vi/ 2.2 and A are small enough. 
Moreover, if [0, T*) denotes the maximal interval of continuation ofu and ifT* < oo then 
ll'a(i) II ^ 4 ^ 2,2 — y oo as t — y T*. 

An analogous result holds for the Navier problem. Note that for this last state¬ 
ment, unlike in the previous results concerning the stationary problems, we have 
allowed the datum h to depend on both space and time coordinates. The proof 
makes use of Banach fixed point theorem and of the potential well techniques and 
can be found in m. 

We have proven more results concerning the parabolic problem (II.Il l but we are 
not going to expose them because of their more technical nature they require more 
preparatory results. They concern the asymptotic properties of the solutions to the 
Dirichlet problem, that converge to zero in the long time provided A = 0 and the 
initial datum is small enough. Also, the blow-up in finite time that takes place for 
both Navier and Dirichlet problems in case the initial datum is large enough. These 
results can be found in IT]. Precisely these results, and in particular those regarding 
the blow-up, motivate in part the study of self-similar solutions to the problems at 
hand. In the following section we outline some facts about this sort of solutions. 

3. The search for self-similar solutions 

In this section we describe an open problem related to the existence of self-similar 
solutions to partial differential equation (fTTI) . For the first time we are going to 
consider this equation set on all the plane instead of on a bounded domain D. Our 
first step is setting A = 0. Next we look for solutions with the following form 
u{r,t) = ^ / {-^), where a and (3 are two real parameters to be fitted in order to 
find a closed ordinary differential equation involving the self-similar variable p = ^ 
and the solution to equation (O expressed as a function of this variable only f{p). 

Substituting our self-similar ansatz in (II.lb we find that this equation adopts the 
closed form 

4/'(i?) - + 4??^/""(b) = 0, 

only if a = 1/4 and /3 = 0. Note the simplicity of this fact possibly makes it one 
of the simplest ways of looking for self-similar solutions. Now we have to provide 
suitable boundary conditions for this ordinary differential equation. Note that by 
its very nature the self-similar variable describes a rotationally invariant solution. 
Therefore we assume the following symmetry conditions on the solution to equa¬ 
tion (EB: dru{0,0,t) = 0 and 9^A^'u(0,0, t) = 0, for a rotationally invariant u, 
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where Ar(-) = ^i9r[r9r(')]- Finally, we impose that both u and Au decay to zero as 
r oo. So this leads to the boundary value problem: 

4/'(r?) - v'^nv) - - 4^V'(^)/"(^) + Sv^riv) + 4r?V""(^) = 0, 

(3.1) /'(O) =/'"(O) = 0, /(r?), /"(r?) ^ 0 as r? ^ (X). 

The existence of solutions to this problem automatically implies the existence of so¬ 
lutions of the form u{x,y,t) = /(p^), to partial differential equation ( II. 11 1 . The 
obvious fact that u = 0 solves (fO with A = 0 and with the assumed boundary 
conditions directly translates into the fact that / = 0 solves boundary value prob¬ 
lem (13.11 1. 

Noting that equation (13.1b does not depend on / but on its derivatives we can ob¬ 
tain a third order ordinary differential equation for g{r]) = The corresponding 

boundary value problem reads 

4g - - 4775' - 477^53' -b Sg '^ g " + ‘ ig ^ g ’" = 0, 

(3.2) 5 ( 0 ) = g"{0) = 0, g'[g) ^0 as 77 ^ 00 . 

Note that this problem is strongly singular and it is to be solved for 77 £ [0,oo). 
Obviously g = 0 solves this boundary value problem. Therefore we will be interested 
in nontrivial solutions. Since one expects the solutions of (13.11) . if any, to behave like 
a Gaussian, one should first try to exclude the cases where f"(0) > 0. We do so in 
the next statement. 

Proposition 3.1. A local solution of the problem 

(3.3) 4g - g'^g - 4gg' - 4g^gg' + Bg^g” + ig^g”' = 0, 

g(0) = 5"(0) = 0, 5'(0) > 0, 

blows up in finite time. More precisely, there exists g > 0 such that g'{g) > 0 forg G (0, 77 ) 
and liiTL^-^^ ( 7 ( 77 ) = -boo. 

Proof. Let g be a local solution of (13.3b . For contradiction, assume that there exists a first 
770 > 0 such that g'{go) = 0. Then g{g) > 0 and g'{g) > 0 for all 77 £ (0, gf). In turn, if we 
rewrite the equation appearing in (13.3b as 

4 ( 77 ^ 5 " - v^g' + vg)' = g^g{^g' + 

this readily shows that {g^g" — g^g' + gg)' > 0 in (0, 770 ]. Since the term inside the bracket 
vanishes at 77 = 0 this yields g^g"{g) — g^g'{g) + 775 ( 77 ) > 0 in ( 0 , 770 ], which we may also 
rewrite as 

//, ^ g'iv) , giv) ^ n f ^ 7 n 1 

g (d)-f >0 for 77 e ( 0 , 770 . 

7] T]^ 
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This proves that {g' — ^)' > 0 and since the term inside this bracket vanishes as p —0, this 
also gives g'{rio) > > 0’ contradicting the characterization of po- We have so proved 

that 


(3.4) 


g'{r]) > 0 for all rj € [0, rj) 


where g is a local solution to (13.3b and -q is the endpoint of its interval of continuation, 
possibly inhnite. Moreover, what we have seen also proves that 


(3.5) 


lira 5 ( 77 ) = +CXD 


in both the cases q < 00 (blow up in hnite time) and 77 = +c» (global solution). 

Now put 77 = e’', g{q) = h{\ogq), h{r) = g{e^). Then the equation in (13.3b reads 

(3.6) 4 (h"'(r) — h"{r) — h'{r)) = 4e^h'(r)/i(r) + (e^’" — 4)/i(r), r G (— 00 , + 00 ), 
while conditions (13.4b and (13.5b become 

(3.7) h'{r)>0 for all r < r, lirn/i(r) =+ 00 , r = logq. 

r—¥r 

For contradiction, assume that r = +c». From (13.6b and (13.7b we infer that 

(h"{r) — h'(r) — h{r))' > e^h'{r)h{r) > 2h'{r)h(r) for all r > log 2 . 

By integrating this inequality over [log 2, r] we get 

h"{r) — h'{r) — h(r) > h{r)'^ + 7 for all r > log 2 

where 7 = /i"(log2) — /i'(log2) — h(log2) — /i(log2)^. By (13.7b we may multiply this 
inequality by h'ir) and maintain its sense: 

h"{r)h'(r) > +h{r)h'{r) + h{r)‘^h'{r)+qh'{r) > h{r)h'{r) + h{r)‘^h'{r)+qh'{r) 

for all r > log 2. Let us now integrate this inequality over [log 2, 7 ’]; we obtain 

s 

where S = i/i'(log 2)^ — ^h{\og 2)^ — 4/i(log 2)^ — 7 /i(log 2 ). By (13.7b we know that there 
exists R > log 2 such that the latter inequality yields h'(r)'^ > jh{r)^ for all r > i?. By 
taking the square root of this inequality we obtain > 5 which, upon integration over 

[i?, r] yields 


1 


> 


r — R 


1 


/7(i?)i/2 - 4 7i(r)i/2 

By letting r +00 we reach a contradiction. 


for all r > R. 


□ 
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Figure 1 . Numerical solutions of the ordinary differential equation cor¬ 
responding to boundary value problem (13.2b . The boundary conditions 
3 ( 0 ) = g"(0) = 0 were employed in this numerical integration. The 
third condition is arbitrary and chosen to be; ^'(O) = —1 (red line), 
( 7 '( 0 ) = —10 (green line) and ^'(O) = —10^ (blue line) in hgure [T(a)l 
g'(0) = —1 (red line), p'(0) = —10^ (green line) and ^'(O) = —10^ (blue 
line) in hgure pXb)l g'(0) = —10^ (yellow line), ( 7 '( 0 ) = —10^ (red line), 
( 7 '( 0 ) = —10® (green line) and ( 7 '( 0 ) = —10^ (blue line) in hgure [T(^ 
g'(0) = -10^ (yellow line), p'(0) = -10® (red line), 5 '(0) = -8.5 x 10® 
(green line) and g'(0) = —2.4 x 10® (blue line) in hgure p^d)] 


In the case g'(0) < 0, we have numerically integrated this differential equation, 
for an example see hgure [TJ Based on Proposition 13.11 and on these preliminary 
numerical experiments we conjecture the following result: 

Conjecture 3.1. There exist no nontrivial solutions to boundary value problem (13.2b . 

It is interesting to highlight one property of the solutions to (13.3b for any g'{0) 
used in the numerical integrations shown in hgure[T] If we assume that g G (^^^[0, 00 ), 
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we have 


9{v) = 9'{0)V + 5"'(0)y + o{v^), g\v) = 9'(0) + g'"(0)y + o(g^), 

9"(9) = 9'"(0)9 + 0 (g), g"'(g) = g'"(0) + o(l) as g ^ 0. 

Substituting this into the equation and simplifying we arrive at 


^5"'(0)-5'(0)^ g^ + o(g^) = 0, 


that implies g'"(0) = 3g'(0)^/8 > 0. Consequently g"(g) > 0 in a right neighborhood 
of 77 = 0 and thus g(g) is convex in a right neighborhood of ry = 0. Note that this 


feature can be observed in the numerical solutions in figure 1(a) 


Let us now make some final remarks. Note first that a proof of coniecture l3.1l would 
not imply the nonexistence of (any kind of) self-similar solutions to partial differential 
equation (fO . The existence of these is actually a question that remains open. Also, 
the trivial solution g = 0 corresponds to a constant /, but the only constant / that 
solves boundary value problem (13.1b is / = 0 due to the extra boundary condition. 
As we have already mentioned, this corresponds to the solution m = 0 of the partial 
differential equation. In connection to this, any nontrivial solution to problem (13.2b 
gives rise to a nontrivial solution / to dO i only if it satisfies an extra cancelation 
condition. 
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